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Abstract 

Motivated by experiments on (Sr,Ca,La)i4Cu2404i, a two-dimensional Ising model with mobile defects and a two- 
dimensional anisotropic Heisenberg antiferromagnet have been proposed and studied recently. We extend previous 
investigations by analysing phase diagrams of both models in external fields using mainly Monte Carlo techniques. 
In the Ising case, the phase transition is due to the thermal instability of defect stripes, in the Heisenberg case 
additional spin-flop structures play an essential role. 
© 2008 Elsevier B.V. All rights reserved. 
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1. Introduction 

The compounds (Sr,La,Ca)i4Cu2404i display inter- 
esting low-dimensional magnetic features arising from 
CU2O3 two-leg ladders and CUO2 chains. Here, we shall 
consider two types of models which have been moti- 
vated by experimental observations on the CUO2 chains 
[1,2,3,4]. 

On the one hand, a simple two-dimensional Ising 
model with mobile defects has been introduced [5] . The 
spins correspond to the magnetic Cu^"*" ions, and the 
defects to those Cu ions which are believed to be spin- 
less due to holes (Zhang-Rice singlets). The defects 
have been shown to form, at low temperatures, nearly 
straight stripes, perpendicular to the CUO2 chains. The 
coherency of the stripes gets lost at a phase transition 
of first order [5,6]. 

On the other hand, Matsuda ei al. [3] proposed a 
two-dimensional anisotropic Heisenberg model, with 
an easy spin axis, which reproduced nicely the mea- 
sured spin-wave dispersions in La5CagCu24 04i. In 
subsequent analyses of this model [7,8], the effect of 
external fields parallel (leading to a spin-fiop phase) 
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and perpendicular to the easy axis on various ther- 
mal quantities has been analysed and compared to 
experimental findings [2,4]. 

In this contribution, we shall extend the previous 
work on both models in external fields. For the Ising 
model with a local pinning potential, the phase dia- 
gram in the (temperature, field)-plane is determined 
at various pinning strengths. In the Heisenberg case, 
the boundary of the antiferromagnetic phase and the 
transition to the spin-flop phase are studied in detail, 
both for the model of Matsuda et al. as well as for a sim- 
pler, more common variant of the anisotropic Heisen- 
berg antiferromagnet [9]. 



2. Ising model with mobile defects 

The Ising model with mobile defects is defined on a 
square or rectangular lattice with one axis defining the 
chain direction (horizontal direction in Fig. 1) [5,6]. 
Each lattice site is occupied either by a spin. Si — ±1, 
or by a defect, Si = 0. Usually, the concentration of 
defects is fixed to be ten percent of the lattice sites, 
as it seems to be the case in La5Ca9Cu24 04i. Defects 
are assumed to be mobile along the chains, keeping a 
minimal distance of two lattice spacings. Neighbouring 
spins are coupled ferromagnetically, J > 0, along the 



0304-8853/08/$ - see frontmatter © 2008 Elsevier B.V. All rights reserved. 



M. Holtschneider et al. / Journal of Magnetism and Magnetic Materials (2008) 1-0 



Fig. 1. Typical equilibrium configuration with zig-zag 
structures in the minimal Ising model, at fcsT/|Ja| = 0.2 
and Ep/\Ja\ =0.1. Only a part of the system with totally 
80 X 80 sites is shown. 
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Fig. 2. Phase diagram of the minimal Ising model 
various pinning strengths, Ep. 
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chains, and antiferromagnetically, Ja < perpendic- 
ular to them. In addition, next-nearest neighbouring 
spins in the same chain separated by a defect interact 
antiferromagnetically, Jo < 0, as suggested by exper- 
iments. Assuming J and Jo to be large compared to 
|Ja|, one arrives at a 'minimal version' of the model, 
where spins along a chain have the same sign between 
two defects, reversing sign at a defect. The only rele- 
vant energy parameter is Ja [5,6]. 

The model, in its minimal and full variants, is known 
to form straight defect stripes, perpendicular to the 
chains, in the ground state, T = 0, with arbitrary sep- 
aration between the stripes. As temperature T is in- 
creased the stripes will meander, tending to keep, on 
average, their largest possible distance due to entropic 
repulsion. Eventually the stripes will break up at a 
phase transition of first order, associated with a pairing 
of defects. We introduce a local pinning with strength 
Ep of the defects at the sites of straight equidistant 
lines perpendicular to the chains. At low temperatures, 
the stripes stay close to the pinning lines and long- 
range antiferromagnetic order is observed. The transi- 
tion remains to be of first order, driven, again, by the 
enhanced pairing of defects [6] . 

Applying now an external field H in the Ising direc- 
tion, the ground state continues to consist of straight 
stripes at sufficiently low fields. At Hi < H < H2, as- 
suming a vanishing or weak pinning potential, the de- 
fects will form zig-zag stripes separating the antiferro- 
magnetic domains, thereby allowing for a non-zero to- 
tal magnetization, see Fig. 1 for a typical equilibrium 
configuration at low temperatures. For H > H2, in 
the minimal model the defects will be paired, with iso- 
lated spins between two neighbouring defects pointing 
opposite to the field direction. In the full model, even 
those spins may be reversed in even stronger fields. 



Upon increasing the temperature, at fixed field H < 
H2, the Ising magnet will undergo a phase transition 
at which the defect stripes, being either straight or of 
zig-zag type, become unstable. Indeed, performing ex- 
tensive Monte Carlo simulations, we determined phase 
diagrams of the Ising model in the (T, J/)-plane, as de- 
picted in Fig. 2 for the minimal variant with different 
pinning strengths, Ep [10]. 

In particular, we find no evidence for a transition as- 
sociated with the zig-zag structure. Albeit there is a 
jump in the total magnetization, and thence a diver- 
gence in the susceptibility, at T = and H = H\ , these 
quantities change smoothly even at very low temper- 
atures. Indeed, small TAg-z&g segments are thermally 
excited well below H\ , with a gradual increase of their 
average length as H is increased [10]. 

A large pinning potential may suppress the zig-zag 
structures in the minimal and full models [10]. 

We conclude that the zig-zag structures are pre- 
sumably not relevant in explaining the transition 
from the antiferromagnetic to the disordered phase 
in (La,Ca)i4Cu2404i [2,4]. However, a melting of ex- 
tended or local defect stripes may play an important 
role in that transition. 



3. Anisotropic Heisenberg antiferromagnet 

Following Matsuda et al. [3] , the magnetic properties 
of La5Ca9Cu24 04i depend on the Cu^"'" ions located 
in the ac-planes, having a centered rectangular geom- 
etry, see Fig. 3. Based on their spin-wave analysis, the 
spins [S — 1/2) of the ions couple along the CUO2 
chains, i.e. along the c axis (vertical direction in Fig. 
3), through nearest neighbour, Jd, and next-nearest 
neighbour, Jc2, exchange constants, with Jd — —0.2 
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Fig. 3. Typical equilibrium configuration of the z-coinpo- 
nents of the spins in the antiferromagnetic phase of a classi- 
cal variant of the anisotropic Heisenberg model of Matsuda 
et. al. [3] without field, at fcsT= 0.5 meV. Only a part of 
the system with totally 40 x 40 sites is shown. 
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meV being antiferromagnetic and Jc2 = 0.18 meV be- 
ing ferromagnetic. Tlie ferromagnetic ordering in tlie 
chains is due to the strong antiferromagnetic inter- 
chain couphngs: Jaci ~ —0.681 meV refers to the two 
nearest neighbours in the adjacent chain, and Jac2 ~ 
—0.3405 meV denotes the couplings to the two next- 
nearest neighbours. 

Importantly, there is an uniaxial anisotropy favour- 
ing alignment of the spins along the b axis. Its contri- 
butions to the different couplings are not known, and 
its total effect may be mimicked in the classical variant 
of the model with spins of length one by a single-ion 
interaction D — —0.211 meV. 

We studied the classical model with external fields 
along the easy axis, H^, and perpendicular to it, Hx, 
doing Monte Carlo simulations. In both cases, one 
encounters an antiferromagnetic phase at low fields, 
see Fig. 3. The complete phase diagrams, as obtained 
from finite-size analyses of the simulational data, are 
shown in Fig. 4. In the case of Hx, the phase tran- 
sition is continuous, being in the Ising universality 
class. More interestingly, in the case of H^, there is 
a spin-fiop-phase with algebraically decaying spin 
correlations. The boundary of the antiferromagnetic 
phase is, at low fields, in the Ising universality class as 
well. As suggested by the behaviour of the Binder cu- 
mulant, the transition between the antiferromagnetic 
and paramagnetic phases eventually becomes of first 
order as the field is increased, with a tricritical point at 
about ksT « 0.79 meV [8]. At even stronger fields, one 
encounters, closeby, at kBT « 0.75 meV [8], a triple 
point between these two and the spin-flop phases. 

Our preliminary Monte Carlo study on the stan- 
dard classical anisotropic Heisenberg model with near- 
est neighbour antiferromagnetic interactions in two di- 
mensions [9] suggests that the phase diagram keeps the 
same topology for that simpler model, as had been pro- 
posed before for its quantum version [11]. 

Note that experiments on La5CagCu24 04i do not 



Fig. 4. Phase diagrams of the anisotropic Heisenberg model 
[3] with fields parallel and perpendicular to the easy axis 
of the spin couplings. 

show a direct transition between the antiferromagnetic 
and spin-flop phases [2,4,7] . This may be due to defects 
which have not been taken into account in the model 
of Matsuda et al. Their possible influence has been 
discussed elsewhere [8]. 
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